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ABSTRACT 

High-quality  measurements  of  the  turbulent  dissipation  rate  e  are  required  to  di¬ 
agnose  field  surfzone  turbulence  budgets.  Quality  control  (QC)  methods  are  pre¬ 
sented  for  estimating  surfzone  e  with  Acoustic  Doppler  Velocimeter  (ADV)  data. 

Bad  ADV  velocity  data  points  are  diagnosed  with  both  the  ADV  signal-strength 
(SS)  and  correlation  (CORR).  The  fraction  of  bad  SS  data  points  (5ss)  depends  in¬ 
versely  upon  the  wave-amplitude  normalized  transducer  distance  below  the  mean 
sea-surface.  The  fraction  of  bad  CORR  data  points  (5coRR  can  be  elevated  when 
(5ss  is  low.  The  (5coRR  depends  inversely  upon  the  wave-amplitude  normalized 
sensing  volume  distance  below  the  mean  sea-surface  and  also  increases  with  in¬ 
creased  wave-breaking,  consistent  with  turbulence  and  bubble  induced  Doppler 
noise.  Velocity  spectra  derived  from  both  “patched”  and  “interpolated”  time  se¬ 
ries  is  used  to  estimate  e.  Two  QC  tests,  based  upon  the  properties  of  a  turbulent 
inertial-subrange,  are  used  to  reject  bad  e  data  runs.  The  first  test  checks  that  the 
vertical  velocity  spectrum’s  power-law  exponent  is  near  —5/3.  The  second  test 
checks  that  a  ratio  R  of  horizontal  and  vertical  velocity  spectra  is  near  one.  Over 
all  (icoRR>  70%  of  patched  and  interpolated  data  runs  pass  these  tests.  However, 
for  larger  (icoRR  >  O-l  (locations  higher  in  the  water  column)  50%  more  patched 
than  interpolated  data  runs  pass  the  QC  tests.  Previous  QC  methods  designed  for 
wave  studies  are  not  appropriate  for  e  QC.  The  results  suggest  that  e  can  be  consis¬ 
tently  estimated  over  the  lower  60%  of  the  water  column  and  >  0.1  m  above  the 
bed  within  a  saturated  surfzone. 
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1.  Introduction 

Surfzone  turbulence  vertically  mixes  momentum,  tracers,  and  sediment.  High  quality  surfzone 
turbulence  measurements  are  critical  to  diagnosing  surfzone  turbulence  energetics.  Measurements 
of  turbulent  dissipation  rate  e,  often  used  to  study  oceanic  turbulence  {e.g.,  Terray  et  al.  1996; 
Gerbi  et  al.  2009),  are  sparse  within  the  surfzone.  Measuring  surfzone  turbulence  is  challenging 
because  breaking  waves  and  strong  currents  exert  powerful  forces  on  instruments,  the  water  and 
sea-bed  both  vary  substantially,  and  high  levels  of  turbulence,  bubbles  and  suspended  sediment 
corrupt  velocity  measurements. 

Acoustic  Doppler  Velocimeters  (ADV)  measures  3-components  of  velocity,  at  high  sampling 
rates,  by  measuring  the  Doppler  shift  of  returned  acoustic  pulses  (Sontek  2004).  ADVs  have  been 
used  to  study  waves  {e.g.,  Thomson  et  al.  2007)  and  mean  circulation  {e.g.,  Apotsos  et  al.  2008) 
in  the  surfzone  and  nearshore.  The  ADV  accurately  measures  Reynolds  stresses  and  turbulent 
velocity  spectra  in  laboratory  flumes  (Voulgaris  and  Trowbridge  2001).  ADVs  have  also  been 
used  to  study  turbulence  in  a  laboratory  surfzone  {e.g.,  Scott  et  al.  2005),  in  estuarine  and  coastal 
{e.g.,  Kim  et  al.  2000)  environments,  and  in  field  surfzone  (Bryan  et  al.  2003)  and  swash-zone 
(Raubenheimber  et  al.  2004). 

The  ADV  sensor  also  returns  the  backscattered  acoustic  signal-strength  (SS)  and  the  corre¬ 
lation  (CORR)  of  successive  pings  {e.g.,  Zedel  et  al.  1996).  Both  SS  and  CORR  are  used  to 
diagnose  ADV  data  quality.  Surfzone  ADV  velocity  measurements  can  be  noisy  with  significant 
amounts  of  bad  data  (Elgar  et  al.  2001).  The  signal  strength  (SS)  depends  upon  the  density  of 
scatterers  {e.g.,  Lohrmann  et  al.  1994).  With  insufficient  scatterers,  SS  is  low  and  the  velocity 
signal  is  unreliable.  Within  the  surfzone,  there  is  generally  no  shortage  of  scatterers  {e.g.,  bubbles 
and  suspended  sediment).  Low  SS  also  occurs  when  the  ADV  sensor  is  exposed  out  of  the  water 
{i.e.,  above  water  level)  or  when  the  scatterer  density  is  so  high  that  the  acoustic  signal  is  absorbed 
or  scattered  {e.g.,  Elgar  et  al.  2005).  Along  beam  correlation  (CORR),  the  coherence  between  the 
Doppler  shift  observed  with  successive  pings,  is  low  (Cabrera  et  al.  1987)  when  scatterers  leave 
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the  sampling  volume  between  pings  or  when  velocity  fluctuates  or  is  sheared  within  the  sample 
volume  (Lhermitte  and  Lemmin  1994).  Low  CORK  lead  to  inaccurate  velocity  estimates  (Zedel 
et  al.  1996).  Low  CORK  also  occurs  in  the  presence  of  significant  number  of  bubbles  (Mori  et  al. 
2007),  and  when  the  ADV  sensing  volume  is  too  close  to  the  bed  (Martin  et  al.  2002;  Elgar  et  al. 
2005),  generally  one  sample  volume  width  above  the  bed  (Raubenheimber  et  al.  2004). 

Elgar  et  al.  (2005)  suggests  flagging  data  points  as  bad  when  the  signal-strength  SS  <  7ss  (SS 
is  an  8-bit  count  0-255)  from  any  of  the  Sontek  Ocean  3  ADV  acoustic  beams.  The  signal-strength 
cutoff  7ss  =  100  was  chosen  by  examination  of  surfzone  data  during  times  when  the  probe  was 
known  to  be  both  in  and  out  of  the  water.  This  733  is  specific  to  the  particular  Sontek  ADV  sensor 
and  other  sensors  may  give  different  733  (B.  Raubenheimer,  personal  communication).  Elgar  et  al. 
(2005)  reject  an  entire  data  run  if  the  fraction  of  bad  SS  data  points  533  >  0.008.  This  empirical 
criteria  is  conservative,  assuring  that  little  bad  data  passes. 

The  ADV  correlation  (CORR)  signal  (ranging  from  0-1.0)  also  is  used  to  diagnose  data  quality 
(Zedel  et  al.  1996;  Sontek  2004).  To  mark  data  points  as  bad,  Elgar  et  al.  (2005)  proposed  a 
correlation  threshold  7corr  of 

7corr  =  0.3  -f  0.4v//.//max  (1) 

where  fs  is  the  sample  frequency  and  /max  is  the  maximum  ADV  sampling  frequency  (/max  = 
25  Hz  for  the  Sontek  Ocean  ADV).  The  upper  (0.7)  and  lower  (0.3)  7corr  limits  are  based  on 
Sontek  (2004)  estimates  for  full  sampling  and  mean  flow,  respectively.  Data  points  where  the 
CORR  <  7coRR  on  any  of  the  three  beams  are  marked  bad.  Unlike  the  ^33  rejection  criteria,  Elgar 
et  al.  (2005)  did  not  propose  a  data-run  rejection  criteria  based  upon  the  fraction  of  bad  CORR 
points  (5corr-  Instead,  two  QC  tests,  based  upon  the  expected  properties  of  the  surface  gravity 
wave  field  in  the  sea-swell  band,  are  used  to  reject  data  runs.  The  first  test  statistic  is  the  pressure  p 
to  (interpolated-)  cross-shore  velocity  (u)  coherence  Cpu,  and  the  second  statistic  is  based  upon  the 
ratio  of  pressure  to  horizontal  velocity  variance  {i.e.,  Guza  and  Thornton  1980).  Thresholds 
for  and  Cpu  tests  were  selected  empirically.  A  large  (5corr  usually  indicated  that  a  data  run 
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would  not  pass  the  and  Cp„  tests,  but  some  eases  with  large  (5corr  (up  to  (5corr  =  0.5)  did 
pass  these  tests  (Elgar  et  al.  2005).  This  ADV  QC  methodology  works  well  for  wave  and  mean 
eurrent  studies  {i.e.,  at  frequeneies  <  0.3  Hz),  but  is  eonstrained  by  the  requirement  of  a  eo-loeated 
and  synehronized  pressure  measurement. 

Turbulent  dissipation  rate  e  estimates  depend  erueially  upon  the  high  frequeney  (>  1  Hz) 
eomponent  of  the  veloeity  speetrum  rather  than  on  sea-swell  band  frequeneies  (~  0.1  Hz).  In 
nearshore  and  surfzone  field  studies  of  e,  ADV  QC  methods  vary.  In  4.5  m  mean  water  depth,  a 
loeation  only  oeeasionally  within  the  surfzone,  ADV  measurements  1  m  above  the  bed  were  used 
to  estimate  e  (Trowbridge  and  Elgar  2001).  ADV  data  quality  eontrol  used  the  manufaeturer’s 
suggested  7corr;  <5corr  levels,  and  bad  data  interpolation  were  not  diseussed.  Instead,  assuming 
unidireetional  and  pure  shallow-water  wave-orbital  motions  with  a  steady  eurrent,  Trowbridge 
and  Elgar  (2001)  showed  that  within  an  inertial  subrange  (at  high  frequeneies  with  no  instrument 
noise) 

(12/21)(P„„(/)  +  P„„(/))  _ 

Pu.M)  ’ 

where  Puu,  Pw,  and  P^^  are  the  eross-shore,  alongshore  and  vertieal  veloeity  speetra,  respeetively 
and  /  is  frequeney.  The  quality  of  e  estimates  was  ensured  by  eheeking  that  a  ratio  R,  based  upon 
the  left-hand-side  of  (2),  was  near  one.  Although  the  assumptions  used  to  derive  (2)  are  not 
generally  valid  within  a  natural  surfzone,  on  average  R  ^  0.8  was  observed,  suggesting  that  a 
turbulent  inertial  subrange  was  present  and  the  resulting  10  m^  s  ^  (Trowbridge  and  Elgar 
2001). 

In  a  study  of  the  e  vertieal  strueture  seaward  of  the  surfzone  (no  depth-limited  wave  breaking) 
in  3.5  m  mean  water  depth  (Eeddersen  et  al.  2007),  the  SS  and  CORR  (1)  thresholds  (Elgar  et  al. 
2005)  were  applied  to  data  from  a  vertieal  staek  of  3  ADVs.  No  eo-loeated  pressure  sensor  was 
present  and  thus  the  Z'^  and  Cpu  tests  (Elgar  et  al.  2005)  eould  not  be  applied.  Eor  all  ADV  data 
runs,  the  maximum  (5corr  =  0.026,  and  ^corr  was  typieally  mueh  less.  Data  flagged  as  bad 
was  interpolated  following  Elgar  et  al.  (2005),  and  the  resulting  R  varied  between  0.8-1. 5,  and 


6 


JOURNAL  OF  ATMOSPHERIC  AND  OCEANIC  TECH 


Volume 


e  varied  between  10“^-3  x  10“"^  s“^.  However,  these  observations  did  not  reach  far  up  in  the 

water  column  as  z' /Hs\g  >  1  {z'  is  the  ADV  distance  below  the  mean  sea-surface  and  ifsig  is  the 
significant  wave  height). 

In  an  study  of  turbulent  energetics  with  white-capping  wave  breaking  in  16  m  water  depth 
(Gerbi  et  al.  2009),  data  runs  with  large  vertical  velocities  or  the  ADV  sensor  too  close  to  the 
surface  {i.e.,  z'/Hsig  <  2)  were  rejected,  and  the  resulting  e  <  10“^  m^  s“^.  In  a  shallow  estuary 
(1.5-3. 5  m  depth)  study  (Jones  and  Monismith  2008)  of  the  e  vertical  structure  with  whitecapping 
wind- waves  (Hgig  between  0. 1-0.6  m),  the  ADV  velocity  QC  methods  were  not  specified.  Esti¬ 
mates  of  e  were  rejected  if  the  vertical  velocity  spectrum  was  not  consistent  with  a  —5/3  power 
law  over  some  frequency  range.  Measurements  were  reported  relatively  high  up  in  the  water  col¬ 
umn  with  z' /Hsig  as  small  as  0.3  and  zjh  (where  2;  is  height  above  the  bed  and  h  is  the  mean  water 
depth)  as  large  as  0.9.  The  resulting  e  was  generally  <  10“"^  m^  but  occasionally  was  as  large 
as  10“^  m^  high  in  the  water  column. 

In  a  study  of  surfzone  e  in  <  3  m  depths  and  with  incident  H^ig  <  0.6  m  (Bryan  et  al.  2003), 
Sontek  Ocean  ADV  data  points  with  SS  and  CORK  below  the  733  =  77  counts  and  tcorr  =  0.7 
were  marked  as  bad.  Data  runs  were  rejected  if  the  fraction  of  total  bad  data  points  >0.1,  resulting 
in  62  of  194  data  runs  being  discarded.  The  data  interpolation  method  was  not  specified.  Data 
runs  were  additionally  rejected  if  the  best-fit  velocity  spectra  power-law  was  not  near  —5/3.  The 
resulting  e  varied  between  10“^-10“^m^  s“^.  Some  retained  data  runs  were  relatively  high  up  in 
the  water  column,  at  time  exceeding  zjh  >  0.7. 

In  a  swash-zone  turbulence  study  (Raubenheimber  et  al.  2004),  two  vertical  stacks  of  (two- 
velocity  component)  ADVs  were  deployed  in  5  cm  and  25  cm  mean  water  depth.  At  the  25  cm 
location,  the  ADV  was  considered  submerged  (from  the  SS  signal)  most  (98%)  of  the  time.  At 
the  5  cm  location,  only  25%  of  the  data  runs  were  considered  as  submerged.  Bad  data  points 
were  removed  following  Elgar  et  al.  (2005).  Data  runs  considered  submerged  rarely  had  (5corr  > 
0.03.  Velocity  spectra  were  calculated  from  the  Eourier  transform  of  the  velocity  auto-covariance. 
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precluding  the  need  for  data-gap  interpolation,  but  perhaps  biasing  the  spectra  estimates.  The 
turbulent  dissipation  rate  e  was  estimated  from  the  high-frequency  spectra  following  Trowbridge 
and  Elgar  (2001).  The  observed  e  were  the  largest  oceanic  e  ever  reported,  up  to  10“^  m^s“^, 
and  were  an  order  of  magnitude  larger  than  the  combined  shear  production  and  depth-normalized 
breaking  wave  energy  flux  gradient  (Raubenheimber  et  al.  2004). 

Here,  surfzone  and  nearshore  ADV  data  are  used  to  examine  and  develop  quality  control 
methodology  for  estimating  surfzone  e.  This  process  also  should  be  applicable  to  open-ocean 
air-sea  boundary  e  studies.  The  surfzone  and  nearshore  field  ADV  observations  from  the  HB06 
field  experiment  (Huntington  Beach  CA,  Fall  2006)  are  described  in  Section  2.  Bad  ADV  ve¬ 
locity  data  points  are  identified  (Elgar  et  al.  2005)  with  both  the  ADV  signal-strength  (SS)  and 
correlation  (CORK)  signals  (Section  3).  The  fraction  of  bad  SS  data  points  (5ss  is  a  function  of 
the  (wave-amplitude)  normalized  ADV  transducer  depth  below  the  surface.  The  bad-SS  data  gap 
statistics  are  used  to  assist  in  identifying  a  (iss  cutoff  to  reject  bad-SS  data  runs.  The  fraction  of 
bad  CORR  points  (5corr  can  be  large  even  with  small  (5ss-  The  resulting  (5corr  is  related  to  both 
the  sensing  volume  distance  below  the  surface  and  the  wave  energy-flux  gradient,  consistent  with 
turbulence-  and  bubble-induced  Doppler  noise  within  the  sensing  volume. 

The  method  for  estimating  e,  the  QC  tests,  and  their  application  are  described  in  Section  4.  The 
two  QC  tests  are  based  upon  the  properties  of  the  turbulent  inertial-subrange.  Velocity  spectra  are 
calculated  from  “patched”  and  “interpolated”  time-series.  At  smaller  (icoRR  (<  OT),  patching  and 
interpolation  give  similar  results.  At  higher  (5corr>  patched  data  runs  are  more  often  consistent 
with  an  inertial  subrange,  and  some  data  runs  pass  with  (5corr  as  high  as  0.4.  The  implications 
of  the  e  QC  method  are  discussed  in  Section  5.  The  interpolated  e  estimates  are  biased  low 
relative  to  patched  e  estimates.  Previous  QC  methods  designed  for  wave  studies  are  shown  to 
be  inappropriate  for  e  QC.  Surfzone  e  estimates  can  be  consistently  made  at  about  1.5  x  wave 
amplitude  below  the  mean  sea-surface,  corresponding  to  the  lower  60%  of  the  water  column  in  a 
saturated  surfzone.  The  results  are  summarized  in  Section  6. 
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2.  The  HB06  Surfzone  ADV  Observations 


a.  HB06  Instrumentation  and  Processing 

Fig. 1 1 

Surfzone  field  observations  were  collected  during  Fall  2006  at  Huntington  Beach  CA  State 
Park  (lat:  33.636  N,  Ion:  -117.969  E)  as  part  of  the  HB06  experiment  (Spydell  et  al.  2009;  Clark 
et  al.  2010;  Omand  et  al.  2010).  A  cross-shore  transect  of  six  instrumented  frames  was  deployed 
spanning  160  m  from  near  the  shoreline  out  to  4  m  mean  water  depth  (Fig.  1).  An  additional 
deployed  instrumented  frame  (between  instruments  1-2)  was  often  buried  and  observations  from 
it  are  not  included  here.  At  each  instrument  location,  the  vertical  coordinate  is  positive  upward 
with  2;  =  0  m  at  the  bed.  The  cross-shore  coordinate  x  is  positive  offshore.  The  instrument  frames 
were  leveled  with  possible  orientation  errors  of  ±3°.  The  tide  range  was  approximately  ±1  m. 
Data  was  collected  for  800  hours  from  14  September  to  17  October  2006. 

Each  instrumented  frame  had  a  buried  pressure  (p)  sensor  and  mounted  downward  looking 
5  MHz  Sontek  Ocean  Probe  ADV  (Sontek  2004)  with  synchronized  data  collection  sampled  at 
8  Hz.  Vertical  instrument  locations  were  GPS  measured  to  within  a  few  cm  relative  to  mean  sea 
level.  The  ADV  measures  3  components  of  velocity  (u,  v,  and  w)  aligned  with  the  coordinate 
system.  The  velocity  range  was  set  to  ±5  ms“\  and  velocities  beyond  this  range  (i.e.,  phase 
wrapping)  were  not  observed.  In  addition,  the  ADV  returns  signal-strength  (SS)  and  correlation 
(CORK)  on  each  of  the  3  beams.  Both  SS  and  CORK  are  given  as  an  unsigned  byte  (0-255  counts) 
and  CORR  is  normalized  to  between  0-1.0.  In  each  hourly  data  run,  the  ADV  sampled  24,578  data 
points  (51.2  min  or  3072  s)  and  subsequently  went  into  bottom  finding  mode  for  the  remainder  of 
the  hour  to  estimate  ADV  transducer  height  above  the  sea-bed  (ztr)  and  bed  location  (relative  to 
mean  sea  level). 


Prom  each  pressure  sensor,  the  mean  sea  surface  location,  mean  water  depth  h,  and  sea- 
surface  elevation  spectra  iPr^n)  were  estimated  hourly.  These  calculations  are  independent  of 
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the  co-located  ADV  velocity  data.  Pressure  spectra  from  buried  sensors  were  adjusted  follow¬ 
ing  Raubenheimber  et  al.  (1998).  From  the  spectra,  significant  wave  height  Hgig  is  calculated 
over  both  sea-swell  (0.03-0.3  Hz,  infragravity  (0.003-0.03  Hz,  frequency  bands 

The  total  (sea-swell  and  infragravity  bands)  significant  wave  amplitude  Osig  is  given  by  Usig  = 

((i/gig)^  -f  (iFg?g)^)^/^/2.  During  the  experiment,  the  incident  varied  between  0.5-1.4  m. 

The  (downward  looking)  ADV  sensing  volume  vertical  location  Za_dv  is  0.18  m  below  the 
transducer  location  Ztj-  {i.e.,  Ztr  =  z^dv  +  0.18  m).  The  ADV  sensing  volume  is  a  approximately 
(0.01  m)^)  cylinder  (Sontek  2004).  During  the  deployment  the  sea-bed  eroded  and  accreted  and 
the  ADVs  were  occasionally  raised  or  lowered  on  the  frames.  At  instruments  1-3,  Zs,dv  varied 
between  0-0.4  m,  and  at  instruments  4-6,  Zg^v  varied  between  0.5-0. 8  m.  Data  runs  with  sensing 
volume  too  close  to  the  bed  (zadv  <  0.03  m)  are  rejected.  The  distance  below  the  mean  sea-surface 
of  the  sensing  volume  and  transducer  (z^j.)  is  given  by  =  h  —  ;2adv  and  z[^  =  h  —  Ztr, 
respectively.  Both  z'^^^  and  zj-j.  are  relevant  because  when  the  transducer  of  a  downward  looking 
ADV  is  exposed  out  of  the  water  the  acoustic  path  is  blocked  even  if  the  sensing  volume  remains 
submerged.  For  an  upward  looking  ADV  this  is  not  a  concern  as  the  sensing  volume  location 
would  be  exposed  first.  For  a  horizontally  mounted  ADV  ztr  =  Zadv 

b.  Example  of  ADV  Data 

Fig.  2 

The  challenges  in  using  surfzone  ADV  data  to  estimate  (high  frequency)  turbulence  parameters 
are  illustrated  with  a  short  (160  s)  time  series  of  ADV  data  (Fig.  2).  In  general,  the  vertical 
velocities  are  small  (\w\  <0.1  ms“^.  Fig.  2a)  as  expected  for  shallow  water  surface  gravity  waves. 

The  signal  strength  is  typically  SS  >  180  counts  (Fig.  2b),  well  above  the  suggested  733  =  100 
counts  cutoff  (Elgar  et  al.  2005).  In  addition,  correlations  generally  are  high  >  0.8  (Fig.  2c),  above 
the  7CORR  =  0.526  (Eq.  1  with  fs  =  8  Hz)  cutoff  (Elgar  et  al.  2001).  However,  occasionally  the 
vertical  velocities  are  quite  large  (exceeding  |t(;|  >  0.8  ms“^).  and  noisy  (e.g.,  near  35  s  in  Eig.  2a) 
when  SS  (Eig.  2b)  or  CORK  (Eig.  2c)  are  low,  falling  below  the  suggested  733  and  7corr  cutoffs. 
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The  SS  has  a  minima  at  an  apparent  noise  floor  of  42  counts.  CORK  can  fall  below  Tcorr  =  0.526 
when  SS  does  not  (i.e.,  near  20  s  in  Fig.  2).  For  the  entire  data  run,  (5ss  =  0.016,  exceeding  the 
Elgar  et  al.  (2005)  (5ss  =  0.008  cutoff,  which  would  result  in  rejection  of  this  data  run.  The  fraction 
of  combined  bad  SS  and  CORR  data  points  (5corr  =  0.045,  exceeding  that  typically  observed  by 
Raubenheimber  et  al.  (2004)  and  Feddersen  et  al.  (2007).  It  is  not  known  whether  this  level  of 
(icoRR  can  be  tolerated  in  estimating  e. 

3.  Quality  Control  (QC)  of  ADV  Data 


a.  Signal  Strength  (SS)  QC  ofSurfzone  ADV  Data 

Fig.  3 

Within  a  data  run,  ADV  data  is  marked  bad  when  the  returned  signal- strength  SS  <  733  at 
any  of  the  3  acoustic  beams  with  733  =  100  counts  (Elgar  et  al.  2005).  With  a  ~  42  count  ADV 
noise  floor  (see  Eig.  2b)  and  a  0.43  dB/count  conversion,  this  733  =  100  count  cutoff  corresponds 
to  a  25  dB  cutoff,  which  is  more  conservative  than  the  15  dB  Sontek  (2004)  recommendation. 
However,  the  resulting  ^33  is  insensitive  to  733  within  the  range  of  80-130  counts  (16-38  dB). 

The  fraction  of  bad  SS  data  runs  (5s3  is  calculated  for  all  the  data  runs.  At  all  instruments,  533 
did  not  systematically  depend  upon  instrument  height  above  the  bed  indicating  that  high  levels  of 
near-bed  suspended  sediment  (Beach  and  Sternberg  1996,  e.g.,  )  does  not  adversely  impact  ADV 
signal-strength.  In  the  nearshore  and  surfzone,  the  sea-surface  fluctuates  due  to  infragravity  and 
sea-swell  surface  gravity  waves  which  can  expose  out  of  the  water  an  instrument  deployed  below 
the  mean  surface.  The  amount  the  ADV  transducer  is  exposed  out  of  the  water  and  thus  533  is 
expected  to  increase  with  smaller  z[^  (the  distance  of  the  downward-facing  ADV  transducer  below 
the  mean  sea  surface)  and  increase  with  larger  significant  wave  amplitude  Ogig.  Reflecting  this, 
^33  is  inversely  related  to  the  normalized  ADV  transducer  depth  z[^l Ogig  (Eig.  3)  with  a  consistent 
relationship  that  collapses  at  all  surfzone  instrument  locations  (1-4).  At  4/agig  =  0.5,  (^33  varies 
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between  0. 1-0.2,  and  for  larger  z'^Ja^ig  >  1  (conceptually,  the  ADV  transducer  below  the  signif¬ 
icant  trough  level),  (5ss  is  much  reduced,  generally  <  0.02.  If  Osig  does  not  include  infragravity 
fluctuations,  the  relationship  between  (5ss  and  z[^/asig  does  not  collapse  as  well,  particularly  near 
the  shoreline  (instrument  1)  where  infragravity  energy  can  be  significant  {e.g.,  Guza  and  Thornton 
1985).  Not  all  large  (iss  result  from  ADV  exposed  out  of  the  water.  Occasionally  (5ss  =  0.1  when 
^(r/®sig  ~  2  (Fig.  3)  and  at  times  (<  2%  of  data  runs)  (5ss  >  0-01  at  instruments  5  and  6,  which 
are  always  well  below  the  mean  surface  {z[^/as\g  >  1.75,2.75,  respectively;  not  shown).  Other 
mechanisms  (lack  of  sufficient  scatterers  or  acoustic  absorption/scattering)  induce  these  moderate 
(5ss-  The  Elgar  et  al.  (2005)  (5ss  <  0.008  criteria  (horizontal  dashed  red  line  in  Fig  3)  rejects  all 
data  runs  with  z[^/as\g  <  1  which  may  be  of  particular  interest  for  turbulence  studies. 

The  functional  form  for  the  ()ss-2:(r/®sig  relationship-boundary  is  approximately  given  by  a 
quartic  relationship, 

logio(5ss)  =  -0.17(4/a,ig  +  c)"  -  1.64(4/asig  +  cf  -  0.32,  (3) 


where  c  =  0.15  (gray-dashed  curve  in  Fig.  3).  Although  not  aesthetically  pleasing,  the  rela¬ 
tionship  (3)  holds  at  all  surfzone  instruments  regardless  of  whether  in  the  swash-zone  (lower  tide 
at  instrument  1)  or  surfzone.  When  designing  surfzone  ADV  deployments,  (3)  yields  a  (5ss  esti¬ 
mate  for  a  downward  facing  ADV.  For  example,  if  a  (5ss  =  0.1  is  tolerable,  then  measurements 
potentially  can  be  made  as  shallow  as  4r/®sig  =  0.48. 

For  each  data  run,  the  bad  SS  data  gaps  are  binned  into  probability  density  functions  (pdf)  of 
data  gap-lengths  from  1/8-60  s.  The  data  gap-length  statistics  dependence  upon  (5ss  is  used  to  help 
determine  criteria  to  reject  data  runs.  The  pdf  maximum  (the  mode)  is  typically  at  or  near  1/8  s 
(1  sample)  for  all  (5ss  (blue  dots  in  Fig.  4).  The  gap-length  means  and  standard  deviations  (std) 
increase  with  increasing  (5ss  (circles  and  asterisks  in  Fig.  4).  At  all  (5ss,  the  data  gap  length  means 
and  std  are  roughly  equal,  and  together  with  a  1  point  mode,  suggests  approximately  exponentially 
distributed  data  gap-lengths.  The  data  gap-length  statistics  dependence  upon  (5ss  is  independent 
of  ADV  location.  For  (5ss  <  0.1,  both  data  gap-length  mean  and  std  are  typically  <  2  s,  and  are 


Fig.  4 
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linear  with  log^o('^ss)-  For  example,  between  the  order  of  magnitude  change  from  (5ss  =  10“^  to 
(5ss  =  10“^,  the  data  gap-length  mean  and  std  only  increase  from  approximately  1  to  2  seconds. 
For  larger  (5ss  (>  0.1),  the  gap-length  means  and  std  increase  rapidly,  suggesting  a  different  nature 
of  exposure  out  of  the  water,  and  that  such  (iss  levels  are  not  tolerable. 

Bad  SS  data  runs  consistently  out  of  the  water  are  rejected  with  a  (5ss  cutoff  (5gg  {i.e.,  data 
runs  with  (5ss  >  (^ss)-  ^  balance  is  sought  in  selecting  (5gg  to  retain  data  runs  higher  in  the  water 
column  of  interest  for  turbulence  studies.  At  typical  mid-surfzone  locations  (2  and  3),  varying 
(5gg  from  0.008  (e.g.,  Elgar  et  al.  2005)  to  0.1  results  in  1/3- 1/2  more  retained  good-SS  data 
runs  with  smaller  z[^/as\g.  As  the  bad  SS  data-gap  statistics  are  still  small  and  increasing  slowly 
(Fig.  4),  the  SS  cutoff  (5gg  =  0.1  is  chosen  to  retain  more  of  the  surfzone  data  runs  within  the 
range  0.6  <  z[^/asig  <  1.5  (Fig.  3)  that  would  otherwise  be  rejected.  The  impact  of  this  choice  is 
subsequently  discussed. 


Fig.  5 


Fig.  6 


b.  Correlation  ( CORK)  QC  of  Surfzone  ADV  Data 

After  rejecting  bad-SS  ((5ss  >  0.1)  data  runs,  the  correlation  (CORK)  QC  is  applied  to  the 
remaining  data  runs.  Data  points  with  CORK  <  Tcorr  on  any  of  the  3  ADV  beams  are  marked 
as  bad,  where  Tcorr  =  0.562  is  given  by  (1)  with  /^  =  8  Hz  (Elgar  et  al.  2005).  Bad  SS  data 
points  are  also  marked  as  bad  CORK.  The  resulting  fraction  of  total  bad  data  points,  denoted 
^CORR,  can  be  significantly  larger  than  (5ss  (Fig.  5).  Even  for  small  (5ss  (<  10“^),  5corr  can 
approach  one,  reflecting  the  different  processes  leading  to  low  signal-strength  (exposure  out  of  the 
water)  and  low  correlation  (Doppler  noise  or  bubbles).  Instruments  2,  3,  and  4,  had  the  strongest 
levels  of  wave  breaking,  consistently  have  the  largest  values  of  (5corr  relative  to  (5ss  (see  legend 
in  Fig.  5).  An  alternative  velocity  QC  algorithm  (“despiking”.  Goring  and  Nikora  2002),  that  uses 
velocity  signal  properties  together  with  a  minimum  CORR  of  0.3,  gives  a  similar  fraction  bad  data 
points  as  (5corr- 


For  examining  (5corr  dependencies,  the  sensing  volume  vertical  location  (Zg^dv  or  as 
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opposed  Ztr  is  the  appropriate  vertical  location  as  Doppler  noise  within  the  sensing  volume  lead 
to  low  correlations  (e.g.,  Lhermitte  and  Lemmin  1994).  The  (5corr  do  not  depend  systematically 
upon  the  elevation  of  the  sensing  volume  above  the  bed  (^adv)-  Wave-breaking  is  a  source  of 
surfzone  turbulence  (George  et  al.  1994;  Bryan  et  al.  2003;  Feddersen  and  Trowbridge  2005)  and 
bubbles  (e.g.,  Deane  and  Stokes  2002)  to  the  upper  water  column.  Thus,  elevated  (icoRR  are  ex¬ 
pected  higher  up  in  the  water  column  and  under  more  intense  breaking  waves.  The  breaking- wave 
turbulence  and  bubble  input  rate  depends  upon  the  wave  energy  flux  gradient  dF/dx,  where  F  is 
the  cross-shore  wave  energy  flux.  The  (5corr  relationship  to  ^^adv/^sig  and  dF/dx  are  examined. 

Assuming  non-reflective,  normally  incident  waves  and  integrating  over  the  sea-swell  band 
(0.05-0.3  Hz),  the  energy  flux  F  is  estimated  at  each  instrument  location  solely  from  pressure  via 

^0.3  Hz 

F  =  g  P,,(f)c,(f)df,  (4) 

2o.05  Hz 

where  g  is  the  gravitation  constant  and  Cg  is  the  linear-theory  group  velocity.  These  wave  energy 
flux  estimates  (4)  are  largely  consistent  with  estimates  derived  from  combined  Pressure-i-ADV  data 
that  take  into  account  non-normal  wave  incidence  and  reflection  (Sheremet  et  al.  2005).  However, 
the  Pressure-!- ADV  based  F  estimates  are  not  independent  of  ADV  data  quality,  and  thus  are 
not  used.  Wave  energy  flux  gradients  dF/dx  are  estimated  at  instruments  1-5  by  differencing 
F  estimates  from  the  neighboring  onshore  and  offshore  instruments.  At  location  I,  F  =  0  is 
assumed  at  the  shoreline. 

Considering  only  good-SS  data  runs,  the  relationship  of  (5corr  to  2;^dv/®sig  (Fig-  6a)  is  anal¬ 
ogous  to  that  for  6ss  (Fig.  3)  with  (icoRR  increasing  with  smaller  ^^adv/^sig-  In  contrast  to  the 
tighter  (5ss  relationship,  the  (icoRR  range  increases  with  ^^adv/^sig-  As  ^^adv/^sig  — ^  1,  the  data 
cloud  becomes  a  nose  and  (5corr  approaches  one  (Fig.  6b).  For  any  data  runs  with  ^^adv/^sig  <  1, 
<icoRR  >0.7  (not  shown).  Note  that  the  non-dimensional  instrument  depths  ^^adv/^sig  are  larger 
than  sig  (Fig.  3)  and  are  not  directly  comparable. 

At  fixed  ^^adv/nsigj  <5corr  is  generally  larger  with  increasing  dF/dx  (note  the  color  stratifica¬ 
tion  in  Fig.  6),  particularly  for  1  <  ^^adv/^sig  <  2.5  (Fig.  6b).  Non-dimensionalized  surfzone 
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dissipation  observations  (e.g.,  George  et  al.  1994)  and  bubbles  (e.g.,  Garret  et  al.  2000)  decay 
with  depth.  The  elevated  (5corr  closer  to  the  surface  and  with  stronger  wave-breaking  is  consistent 
with  small-scale  turbulent-  or  bubble-induced  Doppler  noise  within  the  sensing  volume.  Measure¬ 
ments  closer  to  the  bed  (e.g.,  ^^adv/^sig  =  3),  even  with  large  dF/dx  (red  points  in  Fig.  6a),  can 
have  low  (5corr  (  <  10“^)  as  small-scale  turbulence  and  bubbles  are  reduced  farther  below  the 
surface.  With  the  factors  that  affect  (5corr  understood,  the  dependence  of  the  e  QC  upon  ()corr  is 
examined  next. 

4.  Quality  Control  of  Turbulent  Dissipation  Rate  e 


a.  Calculation  of  e 

Turbulent  dissipation  rate  e  is  estimated  from  observed  (high)  frequency  vertical  velocity  spec¬ 
trum  with  the  Lumley  and  Terray  (1983)  model  that  converts  a  wavenumber  {k)  spectrum  Pww{k) 
to  to  a  frequency  spectrum  Pww{f)  in  a  mixed  wave  and  mean  current  environment.  Variants  of 
this  method  have  been  used  to  estimate  nearshore  e  (Trowbridge  and  Elgar  2001;  Feddersen  et  al. 
2007).  A  Kolmogoroff  inertial  subrange  velocity  wavenumber  spectra  Pww{k)  rs-/  e^/3/c  3/3  (jjje 
to  homogeneous  isotropic  turbulence  (e.g.,  Batchelor  1953)  is  assumed  present.  At  frequencies 
higher  than  the  wave  frequencies,  e  is  derived  from  the  observable  Pww(f)  through  the  model  form 
(Lumley  and  Terray  1983;  Trowbridge  and  Elgar  2001) 

Pww(f)  =  2^27j-)3/2"^"'"'^-^’  ^u,v,w)j  (5) 

where  a  =  1.5  is  Kolmogoroff ’s  constant,  u  and  „  are  the  mean  and  (wave-dominated) 
variance  of  the  3  velocity  components,  and  is  an  integral  over  3D  wavenumber  space  that 
transforms  the  inertial-subrange  /c“3/3  wavenumber  dependence  to  frequency  (Trowbridge  and 
Elgar  2001;  Feddersen  et  al.  2007).  For  all  non-pathological  velocity  mean  and  variances  ~ 
j-5/3  (Qej-bi  et  al.  2009)  meaning  that,  within  an  inertial-subrange,  Pww(f)  ~  As  noise 
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levels  are  lower  for  the  flow  component  parallel  to  the  ADV  orientation,  the  vertical  (parallel  to 
ADV  body)  velocity  spectrum  Pww{f)  is  used  to  estimate  e.  Given  the  observed  Pww{f)  and 
estimated  M^w,  the  estimated  e(/)  are  calculated  via  (5). 

Once  bad  SS  and  CORK  data  points  are  flagged,  gap-free  time  series  are  generated  in  two  ways 
to  calculate  velocity  spectra.  The  first  is  “interpolation”,  following  Elgar  et  al.  (2005),  resulting 
in  a  time  series  denoted  by  Data  gaps  <  1  s  long  (8  data  points)  are  linearly  interpolated 
from  the  good  data  points  bounding  the  gap.  Bad  data  within  the  longer  gaps  is  averaged  together 
and  the  entire  gap  is  set  to  this  constant  average  value.  The  rationale  is  that  velocity  data  noise  is 
unbiased  (as  long  as  instrument  is  in  the  water)  so  that  averaging  the  gap  results  in  a  more  accurate 
mean  current  over  the  gap  (Elgar  et  al.  2005).  The  interpolation  method  acts  analogously  to  a  low- 
pass  filter  biasing  low  the  high-frequency  spectra.  The  second  method  is  “patching”  {e.g., 
which  combines  linear  interpolation  of  short  data  gaps  (<  0.5  s  or  4  data  points)  and  “patching 
together”  of  longer  data  gaps.  Patching  is  illustrated  with  a  discrete  data  sequence 

R’fc)  R’fc+l)  ^k+2:  ■  ■  ■  1  ^k+m:  ^k+m+1 

with  a  bad  data  gap  of  length  m  from  indices  k  +  1  to  k  +  m.  Patching  cuts  out  the  data  from 
the  gap  and  joins  the  good  ends  so  that  Wk+m+i  — >  Wk+i,  reducing  the  time  series  length  by  the 
total  number  of  bad  data  points.  Patching  has  the  potential  for  creating  large  steps  in  the  resulting 
time  series  where  the  data  gap  ends  are  joined,  which  is  expected  to  enhance  (bias  high) 
the  high  frequency  spectrum.  The  interpolation  of  the  shorter  (and  by  far  most  common)  gaps 
reduces  the  amount  of  time-shifting  which  would  otherwise  redistribute  the  spectrum’s  frequency 
distribution.  Quantities  (i.e.,  spectra,  e)  derived  from  patched  and  interpolated  time  series  are 
denoted  with  “(p)”  and  “(i)”  superscripts,  respectively.  Both  patched  and  interpolated  quantities 
are  denoted  with  superscript  “(p,i)”. 

Velocity  spectra  (Puu'Kf)^  Pvv'Kf)^  PwwKf))  calculated  from  the  patched  and  inter¬ 
polated  time  series  using  70  sec  long  data-segments  (detrended,  Hanning  windowed  with  50% 
overlap)  resulting  in  88  degrees  of  freedom.  At  any  frequency,  the  true  spectrum  is  95%  likely  to 
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Fig.  7 


be  found  within  a  factor  of  [0.76, 1.38]  of  the  observed  spectrum.  Analogously,  (/;  u,  ^ 
is  estimated  (see  Feddersen  et  al.  2007)  with  the  velocity  mean  and  variance  from  the  appropriate 
time  series.  Both  and  are  only  calculated  for  good-SS  (i.e.,  (5ss  <  0.1)  data  runs  and 

that  also  pass  a  very  broad  criteria  that  (5corr  <  0.7. 

Both  patched  (e^'’^(/)  and  interpolated  (e^‘^(/))  dissipation  are  estimated  at  Nf  =  56  frequen¬ 
cies  between  1 .2-2  Hz  via 


e"’"(/)  = 


3/2 


(6) 


This  frequency  range  has  been  used  previously  (e.g.,  Trowbridge  and  Elgar  2001;  Feddersen  et  al. 
2007)  as  little  surface  gravity  wave  variance  is  assumed  present  at  these  frequencies.  Consistent 
with  this  assumption,  a  slope-break  is  often  observed  in  velocity  spectra  (e.g.,  near  /  =  0.5  Hz, 
Smyth  and  Hay  2003).  If  the  model  and  inertial-subrange  wavenumber  spectrum  are  correct,  then 
e  should  be  constant  with  /.  At  higher  frequencies  (>  3  Hz),  generally  has  an  approximately 
constant  noise  floor.  Assuming  no  error  induced  by  u  or  ^  error,  the  P^w  spectra  error 
bars  result  in  the  true  found  within  the  interval  [0.66, 1.61]  of  the  observed  e^'’’‘'(/).  Mean 

(frequency-averaged)  dissipation  rate  for  the  data  run  is  calculated  by  averaging  e^'’’'^(/)  over 
all  frequencies.  Alternative  averaging  methods,  i.e.,  e  =  exp[(log(e(/)))  (Feddersen  et  al.  2007), 
result  in  negligible  difference  (typically  1%,  always  <  5%)  to  standard-averaging.  The  e  standard 
error  is  estimated  from  the  variance  of  e(/),  i.e.,  el  =  Var[e(/)]/(iV/  —  1).  Note  that  the  symbol 
e  is  used  to  represent  standard  errors  where  the  symbol  e  is  used  to  represent  dissipation  rate. 

The  resulting  dissipation  estimates  (and  e^‘^)  range  between  10“®  m^  to  3  x  10“^  m^ 
and  increases  with  (5corr  (Fig-  7).  Although  the  ratio  (varying  between  0.03-0.06)  is 

small,  the  e  standard  error  does  not  indicate  whether  a  data  run  e  is  consistent  with  an  inertial 
subrange.  The  increase  with  (5corr  may  be  natural  due  to  the  more  turbulent  the  surfzone  the 
larger  e  and  also  the  larger  (5corr  •  However,  it  is  not  apriori  clear  if  e  estimates  at  a  particular 
^CORR  are  valid,  and  what  (5corr  level  quality  e  estimates  can  be  obtained.  Additional  quality 
control  tests  are  applied  to  reject  data  runs  inconsistent  with  an  expected  inertial-subrange. 
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b.  Example  of  e{f)  frequency  variability 

^IG.  8 

Examples  of  “interpolated”  velocity  spectra  {e.g.,  Pfl;{f))  and  frequency  dependent  dissipa¬ 
tion  (/))  from  two  data  run  at  instrument  4  (with  the  most  intense  wave  breaking  and  strongest 
currents)  are  shown  in  Fig.  8.  In  the  first  example  (Fig.  8a),  the  (5corr  =  0.028  is  moderate, 

Hsig/h  =  0.45,  indicative  of  the  outer-surfzone,  with  observations  relatively  far  (z'^^/osig  =  3.32) 
from  the  surface.  The  horizontal  velocity  spectra  Pfl  +  Pf^  has  a  surface  gravity  wave  peak  (at 
/  =  0.07  Hz),  which  falls  off  rapidly  at  intermediate-frequencies  0.3  </  <  0.7  Hz,  before 
encountering  a  slope-break  at  /  =  0.8  Hz  (red  curve  in  Fig.  8a).  At  higher  frequencies  (0.8- 
3  Hz),  P^f  follows  a  power  law  with  (between  1.2-2  Hz)  best- fit  exponent  /r^‘^(±e|)^)  =  —1.74 
(±0.12,  where  is  the  standard  error  of  /r^‘^),  close  to  the  theoretical  Kolmogoroff  p  =  —5/3 
inertial-subrange  value  (compare  dashed  green  to  thin  blue  curves  in  Fig.  8a).  Consistent  with 
the  best-fit  near  —5/3,  the  estimated  e^'ff)  are  relatively  constant  in  frequency  (black  curve 
in  Fig.  8c).  The  mean  dissipation  =  1.03  x  10“^  m^  (blue  dashed  line  in  Fig.  8c)  and  the 
best-fit  slope  of  e^‘'(/)  with  /,  —9.0  x  10“®  m^  Hz“^  (dotted  line  in  Fig.  8c),  is  statistically 
indistinguishable  from  zero.  The  (/)  95%  confidence  limits  (shaded  gray  region  in  Fig.  8c)  also 
encompasses  the  more  than  95%  of  the  time.  The  patched  spectrum  P/f/,  power-law  exponent 
^(p)  _  ^  0.11,  e^'^ff)  slope  (=  —9.4  x  10“®  m^  Hz“^),  and  =  1.17  x  10“"^  m^  are 

close  to  the  respective  interpolated  quantities.  That  are  near  —5/3  suggests  the  presence  of 
an  inertial-subrange  and  a  quality  e  estimate. 

The  second  example  has  a  larger  Hsig/h  =  0.55  indicative  of  the  inner-surfzone,  measure¬ 
ments  closer  to  the  surface  (^^adv/^sig  =  2.17)  and  larger  (5corr  =  0.255  (Fig.  8b).  Although  the 
velocity  spectra  is  consistent  with  pressure  over  the  sea-swell  band  (e.g.,  Cpu  =  0.91),  at  higher 
(1-3  Hz)  frequencies,  Pf/  +  Pf/  is  not  monotonic  (red  curve  in  Fig.  8b)  and  the  P/j/,  spectra 
falls  off  too  rapidly  with  frequency  (power-slope  of  =  —2.18  ±  0.09)  for  an  inertial-subrange 
(compare  blue  to  dashed-green  curve  in  Fig.  8b).  The  patched  =  —2.13  ±  0.13  is  simi¬ 
lar  to  Consistent  with  —5/3,  e^‘'(/)  is  reduced  with  with  increasing  frequency.  The 
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e®(/)  best-fit  slope,  —2.1  x  10“^  (dotted  line  in  Fig.  8d),  is  significantly  differ¬ 

ent  from  zero.  Note  that  the  e^‘'(/)  error  bars  (shaded  region  in  Fig.  8d)  always  encompass  the 
g(i)  _  4  23  X  10“^  m^  s“^,  indicating  that  this  type  of  test  to  reject  e  estimates  is  insufficient. 
Furthermore,  the  patched  and  interpolated  quantities  are  not  consistent  as  =  1.82,  all  to¬ 

gether  indicating  that  this  data  run  is  inconsistent  with  a  turbulent  inertial-subrange  and  that  this  e 
estimate  should  be  rejected. 

c.  Application  of  e  QC  Tests 

Two  independent  QC  tests,  based  upon  the  expected  presence  of  an  turbulent  inertial-subrange, 
are  applied  to  the  patched  and  interpolated  data  runs  and  evaluated  as  a  function  of  (5corr-  Data 
runs  that  do  not  pass  both  tests  are  considered  inconsistent  with  a  turbulent  inertial  subrange  and 
their  e  estimates  are  rejected.  First,  the  Pwwif)  power-law  exponent  p  is  tested  for  consistency 
with  —5/3.  Second,  a  ratio  of  horizontal  to  vertical  velocity  spectra  is  required  to  be  near  one. 
These  tests,  examining  the  velocity  spectra  frequency  variation  and  a  bulk  (frequency-integrated) 
quantity,  are  examined  separately. 

1)  Spectra  Power-law  Exponent  Consistent  with  an  Inertial  Subrange 

Fig.  9 

For  each  data  run,  the  (patched  and  interpolated)  best-fit  exponents  (with  error  bars 
±e/’’'>)  are  estimated  by  a  least-squares  fit  of  log(P/f4j')  with  log(/)  over  frequencies  1.2-2  Hz,  as 
in  the  case  examples  (Fig.  8).  fit  standard  error  (e/^^  and  e/^)  typically  vary  between  0.09  and  0.15. 
The  estimated  generally  vary  between  —1  and  —2.4  (dots  in  Fig.  9a),  although  the  range  spans 
[—4,  0].  At  all  (5corr>  the  binned-means  are  close  to  —5/3  (diamonds  in  Fig.  9a),  suggesting 
that  often  an  inertial-subrange  is  present  and  that  the  Lumley  and  Terray  (1983)  model  for  convert¬ 
ing  wavenumber  to  frequency  spectra  often  is  applicable  at  all  ^corr  levels.  The  binned-stds 
are  generally  near  0.35  and  do  not  vary  systematically  with  (5corr  (vertical  lines  in  Fig.  9a).  At 
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<^coRR  <  0.1,  the  interpolated  and  patched  are  similar  (Fig.  9b).  At  larger  (5corr>  the 
binned-means  are  consistently  <  —5/3  and  decrease  with  larger  (5corr-  These  steeper  spectral 
slopes  are  an  artifact  of  the  “interpolation”  scheme  which  at  higher  (icoRR  increasingly  reduces 
high  frequency  energy  (i.e.,  is  a  low-pass  filter). 

The  consistency  of  the  estimated  and  with  —5/3,  as  expected  in  an  inertial-subrange 
in  a  wave-current  environment  (Gerbi  et  al.  2009),  are  tested  to  reject  data  runs.  An  analogous  test 
examines  whether  the  e(/)  best-fit  slope  with  /  is  consistent  with  zero.  Applying  either  test  gives 
similar  results,  and  as  the  (i  test  is  more  familiar  {e.g.,  Bryan  et  al.  2003;  Jones  and  Monismith 
2008)  it  is  applied  here.  As  the  log-spectra  are  not  Gaussian,  the  least-squares  standard  errors 
Efj,  are  approximate,  and  rigorous  statistical  tests  on  and  can  not  be  applied.  Instead,  the 
quasi-heuristic  criteria  is  used  where  a  data  run  is  rejected  if  the  jj,  fit-skill  <  0.5  or  if  the  best-fit 
jj,  fall  outside  of  the  region 

jx  —  2e:^  —  A  <  — 5/3  <  /i  -f  2e:^  -f  A,  (7) 

where  A  =  0.06.  Allowing  non-zero  A  gives  the  test  (7)  leeway  given  the  uncertainty  of  the 
underlying  distribution.  If  the  fx  estimates  were  Gaussian,  then  A  =  0  would  correspond  to  95% 
confidence  limits  and  as  typically  ~  0.12,  A  =  0.06  corresponds  to  99%  confidence  limits. 
In  general  the  /i^P’‘^-fit  skill  was  high.  Only  1.3%  and  0.8%  of  the  patched  and  interpolated  data 
runs,  respectively,  were  rejected  due  to  low  skill.  The  first  case  example  with  =  —1.67  ±  0.13 
(Fig.  8a)  passes  the  test  (7),  whereas  the  second  example  with  fx^'^  =  —2.14  ±  0.09  (Fig.  8b)  fails. 
This  criteria  (7)  is  applied  separately  to  all  and  for  the  good-SS  data  runs. 

The  good-yU^p^  data  runs  (passing  the  test  Eq.  7)  generally  fall  within  the  range  —1.9  <  ix^'^'’  < 
—1.4  (Fig.  9c).  The  good-/r‘^'"^  binned-means  are  very  close  to  —5/3  (diamonds  in  Fig.  9c),  except 
for  (icoRR  >  0.4,  and  the  good-yU^p^  binned-stds  are  reduced  to  around  0.14  (relative  to  0.35  in 
Fig.  9a).  For  (5corr  >  0.2,  the  fit-errors  e^(p)  are  approximately  50%  larger  than  at  smaller 
^coRRj  allowing  larger  yu^p'  deviation  from  —5/3  to  pass  the  test  (7).  The  fx^'^  that  pass  (7)  have  a 
<^coRR  dependence  similar  to  yU^p^  (compare  Fig.  9d  to  9c). 
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Fig.  10 


At  various  (5corr>  between  50-80%  of  the  good-SS  data  runs  pass  the  (red  triangles 

in  Fig.  10).  In  total  For  (5corr  <  0.1,  the  number  of  good-/r  data  runs  is  basically  the  same 
for  patching  and  interpolation,  although  generally  is  slightly  greater  than  iV®  (compare  dia¬ 
monds  to  triangles  in  Fig.  10).  At  ()corr  <  10“^,  80%  of  patched  and  interpolated  data  runs  pass 
the  /i  test.  Over  all  ^corr,  71%  and  68%  of  patched  and  interpolated  good-SS  data  runs  pass  the 
/i  test.  However,  at  larger  (5corr  (>  0.1),  is  more  clearly  reduced  relative  to  (compare 
diamonds  to  triangles  in  Fig.  10),  as  are  biased  low  (i.e.,  Fig.  9b),  and  for  binned  (5corr  >  0.1, 
the  total  is  31%  greater  than  the  total 


2)  Ratio  of  Horizontal  to  Vertical  Velocities  Consistent  with  an  Inertial 
Subrange 

Fig.  II 

Although  the  power-law  exponent  /r  test  (7)  rejects  many  data  runs,  some  data  runs  pass  with 
^CORR  as  high  as  0.6.  To  further  test  the  data  runs,  the  second  e  QC  test  examines  the  rela¬ 
tionship  between  horizontal  and  vertical  velocity  spectra  within  an  inertial  subrange.  Previously 
(Trowbridge  and  Elgar  2001;  Feddersen  et  al.  2007),  the  estimated  e  reliability  was  determined  by 
checking  that  the  ratio  i?  ~  1,  where  R  is  based  upon  (2),  and  is  defined  as 

(12/21){/=/=>(P„„(/)  +  P„„(/)-noise)> 

where  ()  represents  a  frequency-average  between  1.2-2  Hz,  and  “noise”  is  the  Puu  +  Pw  ADV 
noise  level  averaged  between  3.1-4  Hz.  For  all  good-SS  data  runs,  and  are  calculated 
via  (8)  from  the  patched  and  interpolated  velocity  spectra,  respectively.  Although  i?  =  1  is  not 
strictly  required,  as  the  assumptions  that  go  into  (2,8)  are  violated,  the  R  dependence  upon  (5corr 
is  examined,  and  R  limits  are  used  to  reject  data  runs  inconsistent  with  an  inertial-subrange. 

First,  consider  the  good-SS  ((5ss  <  0.1)  patched  and  interpolated  data  runs  (Fig.  lla,b).  For 
^CORR  <  10“^,  both  R^'^^  and  are  generally  near  one  (binned-means  between  0.9-1. 2)  but 
vary  at  fixed  (icoRR  with  occasional  outliers  of  >  3.  For  (5corr  <  10“^,  and 
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binned-means  increase  linearly  with  (5corr>  the  binned-std  also  increase,  and  generally  <  2. 
At  larger  (5corr  (>  0.2),  both  and  are  typically  >  2  and  both  binned-means  and  stds 
increase  rapidly  (Fig.  lla,b).  Considering  the  subset  of  good-//  (that  pass  Eq.  7)  data  runs,  the 
overall  and  R‘^''’  dependence  upon  (5corr  (Fig-  1  lc,d)  is  qualitatively  similar  that  for  the  good- 
SS  data  runs  (Fig.  lla,b).  The  good-/i  data  runs  removes  many  of  the  outliers  resulting  in 
binned-means  closer  to  one  and  much  smaller  binned-std.  Thus  the  fj,  and  R  tests  overlap,  as  both 
test  for  an  inertial- subrange.  At  larger  (icoRR  (>  0.1),  the  good-/r  and  binned-means 
increase  more  slowly  than  the  good-SS  R^'^’'\ 

Although  (2,8)  are  not  strictly  valid,  that  and  7?^'^  are  near  one  independently  indicates 
that  an  inertial- subrange  is  often  present.  Based  upon  this,  data  runs  are  rejected  which  do  not 
satisfy  the  heuristicly  chosen  criteria  0.5  <  <  2  (horizontal  dashed  lines  in  Fig.  11).  For 

<^coRR  <  0.1,  very  few  additional  data  runs  are  rejected  with  this  test.  At  (icoRR  >  0.2,  many 
7^(p,i)  >  2,  outside  of  the  selected  heuristic  range.  For  (icoRR  >  0.1,  there  are  1.5  x  the  number  of 
good-/i-7?  patched  versus  interpolated  data  runs  (235  versus  154)  and  the  maximum  (5corr  =  0.58 
and  (icoRR  =  0.34  for  patching  and  interpolation,  respectively.  Thus,  at  higher  (5corr>  patching  is 
more  often  consistent  with  an  inertial-subrange  of  turbulence  and  is  preferable  to  interpolation. 

5.  Discussion 


a.  Ratio 

The  ratio  dependence  upon  (5corr  is  examined  to  determine  their  consistency.  Ideally, 

the  ratio  e'^Ye''^  =  1.  Considering  all  good-SS  data  runs,  the  ratio  depends  upon  (5corr 

(Fig.  12a)  and  generally  >  1,  as  expected.  For  (5corr  <  0.01,  the  are  near  one 

(Fig.  12),  and  the  choice  of  patching  or  interpolation  does  not  impact  e.  For  (icoRR  <  0.1,  the 
binned-mean  slowly  increase  as  does  the  scatter,  but  almost  always  <  1-5-  At 
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Fig.  12 


larger  (5corr  (>  0.2),  the  the  binned-mean  increases  rapidly  as  does  the  scatter.  For  the 

subset  of  good-/i''’’-i?''’’  data  runs,  scatter  is  reduced  (Fig.  12b)  relative  to  only  good- 

SS  data  runs  (Fig.  12a).  However,  the  trend  of  increasing  ^1  larger  (5corr  (>  0.1)  is 

still  present.  For  the  good-/r^'’^-i?‘^'’'  data  runs,  deviations  of  from  1  are  due  to  the  effects 

of  interpolation/averaging  over  longer  gaps  at  higher  (icoRR-  However,  rarely  is  >  2. 

As  e  typically  varies  over  orders  of  magnitude  (here  almost  3  orders  of  magnitude),  the  factor 
of  2  difference  between  U'"’  and  further  indicates  that  the  good-SS-/r^'’^-i?^'’^  patched  are 
accurate. 


b.  Relationship  to  the  p-u  Coherence  Test 

Previous  surfzone  ADV  QC  methodology  (Elgar  et  al.  2001,  2005)  were  designed  for  wave 
and  current  studies  (frequencies  <0.2  Hz),  not  for  estimating  e  (frequencies  between  1-2  Hz).  For 
example,  using  a  synchronized,  co-located  pressure  measurement,  Elgar  et  al.  (2005)  require  that 
the  sea-swell  band  spectral  coherence  (Gpu)  between  p  and  m  is  >  0.9,  based  upon  the  expectation 
that  surfzone  wave  directional  spread  <  25°.  Although,  many  ADV  based  surfzone  and  air-sea 
boundary  e  studies  did  not  have  synchronized  and  co-located  pressure  measurements  (Bryan  et  al. 
2003;  Eeddersen  et  al.  2007;  Jones  and  Monismith  2008;  Gerbi  et  al.  2009),  such  measurements 
were  made  during  HB06  and  the  relationship  between  QC  criteria  and  the  inertial-subrange 
QC  criteria  is  explored. 

Here,  Cp„  is  calculated  at  all  good-SS  data  runs  as  the  sea-swell  band  average  of  the  sea-surface 
elevation  spectrum  weighted  cross-spectral  p-u  coherence  Cp„,  i.e., 

rO.SHz  „  (  f\Af 

Jo.05  Hz 

where  Cpu{f)  is  the  spectral  p-u  coherence  calculated  with  u^'\  and  is  the  (depth-corrected) 
sea-surface  elevation  spectrum.  Eor  the  good-SS  data  runs,  Cpu  varies  between  0. 8-1.0  and  is 
largely  independent  of  (5corr  (Fig-  13a).  The  Cpu  >  0.9  test  (dashed  line  in  Eig.  13a)  is  failed  by 
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33%  of  the  good-SS  data  runs.  The  data  runs  have  a  similar  Cpu  distribution  with 

f^coRR  (Fig-  13b)  to  the  good-SS  data  runs  (Fig.  13a).  The  Cpu  >  0.9  test  is  failed  by  31%  of  these 
good-/i''’'-i?''’'  data  runs.  The  bad-/i^'"^-i?^''^  data  runs  also  have  a  similar  (5corr  dependence  and 
fail  the  Cpu  >  0.9  test  34%  of  the  time  (not  shown).  Thus,  the  Cpu  >  0.9  test  is  equally  likely 
to  pass  or  fail  for  both  good  and  bad  estimates.  This  applies  for  other  Cpu  thresholds  from 
0.8- 1.0,  demonstrating  that  the  Cpu  test  is  not  appropriate  for  quality  controlling  e. 


Fig.  13 


c.  Vertical  Distribution  of  Good  Data  Runs 

Fig.  14 

The  vertical  distribution  of  the  remaining  good  data  runs  is  examined  to  determine  where  in  the 
water  column  e  can  be  estimated.  For  all  good-SS  data  runs,  weakly  decreases  with  smaller 
^ldv/®sig  (gray  dots  in  Fig.  14a).  However,  the  good-yu^p’  data  runs  (red  dots  in  Fig.  14a),  are  in¬ 
dependent  of  ^^adv/^sig’  consistent  with  absence  of  a  good-yU^p^  and  ^corr  relationship  (Fig.  9a,c). 

For  all  good-SS  data  runs,  tends  to  one  at  larger  and  generally  increases  with 

increased  scatter  at  smaller  ^^adv/^sig  (Fig.  14b),  consistent  with  the  relationship  between  and 
^CORR  (Fig.  11).  The  good-yU^p^  values  of  also  follow  this  pattern  with  ^^adv/^sig  (red  dots  in 
Fig.  14b).  Between  1.2  <  z'^-^/osig  <  2,  about  half  of  the  good-yU^'"^  data  runs  are  additionally 
rejected  by  the  limits.  (Fig.  14b).  At  ^^adv/^sig  >  2  where  most  good-yU^p^-i?^''^  data  runs 
are  concentrated,  the  cutoff  (dashed  lines  in  Fig.  14b)  reject  only  a  few  additional  data  runs. 

Although  there  are  a  few  good- if data  runs  as  shallow  as  ^^adv/^sig  =  1-2,  can  only  be 
consistently  estimated  at  ^^adv/^sig  >  When  transformed  from  sensing- volume  (2:adv)  trans¬ 
ducer  coordinates,  this  results  in  a  limit  of  zf/osig  >  1,  largely  corresponding  to  a  stricter  SS 
cutoff  of  Sss  <  10“^  (Fig.  3),  near  the  Elgar  et  al.  (2005)  cutoff  of  Sss  <  0.008. 

The  vertical  velocity  power- law  exponent  yU  has  been  observed  to  transition  from  near  —5/3 
to  —1  within  0.1-0.15  m  above  the  bed  (Smyth  and  Hay  2003)  as  turbulent  eddies  become 
anisotropic.  At  2:adv  <  0.1  m  (sensing  volume  within  0.1  m  of  the  bed),  the  binned-mean  yu^p^ 
deviates  from  —5/3  and  approaches  —1  (not  shown).  In  addition,  the  values  become  larger. 
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consistent  with  anisotropic  eddies.  Few  (22  out  of  146,  15%)  data  runs  passed  both  QC  tests  at 
Za.dv  <  0.1.  At  Zg^dv  >  0.1  m,  no  change  in  or  was  observed.  Thus  z^dv  =  0.1  m  is  a  lower 
near-bed  limit  on  where  can  be  estimated. 

The  ^^adv/^sig  >  1-5  limit  is  useful  in  designing  an  open-ocean  air-sea  boundary  layer  study 
(e.g.,  Gerbi  et  al.  2009).  The  surfzone  is  a  region  of  overlapping  surface  and  bottom  boundary 
layers,  and  from  this  alone  the  water  column  range  where  can  be  estimated  is  not  clear.  Within 
a  saturated  (self-similar)  surfzone  where  Hgig  =  'yh,  where  7  ~  0.5  (Raubenheimber  et  al.  1996), 
the  ^^adv/^sig  ~  1-5  limit  results  in  an  water  column  limit  of  z^dv/h  <  0.6.  Thus,  turbulent 
dissipation  rate  e  can  be  consistently  estimated  in  the  lower  60%  of  the  water  column  and  more 
than  0.1  m  above  the  bed  within  a  saturated  surfzone, 

6.  Summary 

A  quality  control  methodology  for  estimating  surfzone  turbulent  dissipation  rate  e  from  ADV 
observations  is  presented  and  applied  to  HB06  experiment  data.  First,  ADV  velocity  measure¬ 
ments  are  quality  controlled  using  the  ADV  backscattered  signal- strength  (SS)  and  correlation 
signal  (CORK)  to  identify  bad  velocity  data  points.  The  fraction  of  bad  SS  data  points  (5ss  in¬ 
creases  inversely  with  the  (wave-amplitude)  normalized  ADV  transducer  distance  to  the  mean 
sea-surface,  consistent  with  exposure  out  of  the  water  as  the  dominant  reason  for  bad  SS.  Based 
on  statistics  of  the  data-gap  length,  a  liberal  cutoff  criteria  of  (5ss  >  0.1  is  preliminarily  chosen 
to  reject  data  runs.  The  fraction  of  bad  CORR  data  points  (5corr  can  be  significant  even  when 
(5ss  is  small.  The  5corr  is  a  function  of  both  the  (wave-amplitude)  normalized  ADV  sensing  vol¬ 
ume  distance  below  the  mean  sea-surface  and  also  the  wave-energy  flux  gradient,  consistent  with 
turbulence-  and  bubble-induced  Doppler  noise. 

Turbulent  dissipation  rate  e  is  estimated  from  vertical  velocity  spectra  derived  from  both 
patched  and  interpolated  time  series.  Two  QC  tests,  based  upon  the  properties  of  the  expected 
turbulent  inertial-subrange  are  applied  to  reject  bad  e  data  runs.  The  first  test  uses  the  vertical 
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velocity  spectrum’s  power-law  exponent  /i,  expected  to  be  —5/3  in  an  inertial-subrange.  The  sec¬ 
ond  test  checks  that  a  ratio  R  of  horizontal  and  vertical  velocity  spectra  band  is  consistent  with  an 
inertial  subrange.  For  (5corr  <  0.1,  between  60-80%  of  patched  and  interpolated  data  runs  pass 
these  tests.  At  larger  (icoRR  (>  0.1),  50%  more  patched  than  interpolated  data  runs  pass  the  tests, 
and  patched  data  runs  are  used.  Of  the  remaining  data  runs,  the  ratio  of  patched  to  interpolated 
dissipation  is  generally  near  one.  Prior  surfzone  ADV  QC  methodologies  designed  for 

wave  studies  (frequencies  <  0.2  Hz)  have  no  predictive  skill  in  rejecting  bad  e  data  runs.  The 
resulting  good  data  runs  distributed  at  normalized  vertical  locations  z'^^/asig  >  1.5.  This 
suggests  that  turbulent  dissipation  rate  can  be  consistently  estimated  over  the  lower  60%  of  the 
water  column  and  >  0.1  m  above  the  bed  within  a  saturated  (self-similar)  surfzone. 
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Figure  Captions 


Fig.  1.  HB06  cross-shore  depth  transect  versus  distance  from  the  mean  shoreline  (blue  curve). 
The  instrumented  frame  locations  are  given  by  the  circles  and  numbered  1-6.  An  additional 
instrumented  frame,  located  between  1-2,  was  often  buried,  and  is  not  considered  here.  The 
typical  tide  range  is  shown  with  the  horizontal  dashed  red  lines. 


Fig.  2.  Example  ADV  measured  (a)  vertical  velocity  w,  (b)  signal-strength  (SS),  and  (c)  correla¬ 
tion  (CORR)  versus  time.  This  160  s-long  data  segment  is  from  instrument  1  (see  Fig.  1)  on  18 
September  0500.  The  red  dashed  horizontal  line  in  (b)  is  the  suggested  Elgar  et  al.  (2005)  cutoff. 
In  (c)  the  correlation  cutoffs  (  1)  for  /^  =  8  Hz  (tcorr  =  0.526,  red-dashed)  and  the  mean  flow 
/s  =  0  Hz  (tcorr  =  0.3,  magenta)  are  shown.  The  water  depth  h  =  0.57  m,  H^ig  =  0.30  m, 
^^adv  =  0.13  m,  and  =  0.26  m. 

Eig.  3.  Eraction  of  bad  SS  data  points  (5ss  versus  z'^^jasig  at  instruments  1-6  (see  legend)  where  z[^ 
is  the  distance  of  the  ADV  transducer  below  the  mean  sea  surface  and  a^ig  is  the  significant  wave 
amplitude.  Note  that  no  instrument  6  data  points  are  present  in  this  axes  range.  The  horizontal 
dashed-line  is  the  (5ss  =  8  x  10“^  cutoff  for  discarding  a  data  run  (Elgar  et  al.  2005).  The  black- 
dashed  curve  is  the  proposed  scaling  (3)  based  upon  the  data. 


Eig.  4.  Mode,  mean,  and  standard  deviation  of  bad  SS  gap  lengths  versus  ^ss  at  all  instruments. 
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Fig.  5.  (5corr  versus  6ss  at  all  instruments  (see  legend).  The  vertical  dashed  line  indicates  the 
6ss  =  0.1  cutoff. 

Fig.  6.  (a)  The  5corr  versus  ^^adv/^sig  at  instruments  1-5  where  log^o  of  the  energy-flux  gradient 
dF/dx  (m^  s“^)  is  colored.  Only  good  SS  data  runs  that  satisfy  the  (5ss  <0.1  criteria  are  shown. 
Data  runs  with  (5corr  >  0.7  are  not  shown.  Panel  (b)  is  a  blow  up  of  the  nose  regions  in  panel  (a). 


Fig.  7.  e'p'  versus  5corr  for  good-SS  data  runs  at  all  instruments  (see  legend  in  Fig.  5). 

Fig.  8.  Example  of  vertical  velocity  spectra  Pwl,{f)  versus  frequency  /  at  instrument  3  that  are  (a) 
good  (time  =  240,  (5corr  =  0.028,  h  =  1.83  m,  i7sig  =  0.83  m,  ^^adv/^sig  =  3.32,  and  Cpu  =  0.86 
and  (b)  bad  (time  =  280,  ()corr  =  0.255,  h  =  1.03  m,  i7sig  =  0.57  m,  ^^adv/^sig  =  2.17,  and 
Cpu  =  0.91.  The  green  dashed  line  is  a  -5/3  power  slope,  (c)  and  (d)  are  the  versus  /  over  a 
narrower  frequency  range  that  correspond  to  (a)  and  (b).  The  solid  black  line  is  e^‘^(/),  the  shaded 
region  are  the  error  bars  (derived  from  blue  dashed  line  is  and  the  red  dotted  line 

is  the  linear  best  fit  slope. 

Fig.  9.  Vertical  velocity  spectra  Puiw  power-law  exponent  (a)  and  (b)  versus  (5corr  for 
all  good  SS  data  runs  ((5ss  <  0.1).  In  (c)  and  (d)  versus  (5corr  for  cases  that  additionally 
also  pass  the  —5/3-exponent  test  (7).  The  individual  data  points  are  represented  as  gray  dots  and 
the  binned- means  and  stds  are  shown  as  diamonds  and  vertical  bars,  respectively.  The  horizontal 
dashed  black  line  is  the  —5/3  slope.  Note  the  change  in  vertical  scale  between  (a,b)  and  (c,d). 
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Fig.  10.  Binned  number  N  of  good-SS  (N,  circles),  good-/r^'"^  red  triangles),  and  good-/i^‘^ 
green  diamonds)  data  runs  versus  ()corr- 


Fig.  11.  (a)  and  (b)  versus  ()corr  at  all  instruments  for  the  good  SS  ((5ss  <  0.1)  data 
runs,  and  (c)  and  (d)  R^''>  versus  (5corr  for  good-zi^^'  and  good-zr^‘'  data  runs,  respectively. 
Individual  data  points  are  represented  as  gray  dots  and  the  binned-means  and  stds  are  shown  as 
diamonds  and  vertical  bars,  respectively.  The  dash-dot  horizontal  line  represents  R  =  1.  The 
horizontal  dashed  lines  are  the  R  =  2  and  R  =  0.5  cutoffs,  respectively.  Note  the  change  in 
vertical  scale  between  upper  and  lower  panels. 

Fig.  12.  versus  (5corr  at  all  instruments  for  (a)  good  SS  data  runs  ((5ss  <  0.1)  and  (b) 

good-/i^'’'-i?^'’'  data  runs.  The  individual  data  points  are  represented  as  gray  dots  and  the  binned- 
means  and  stds  are  shown  as  diamonds  and  vertical  bars,  respectively. 

Fig.  13.  Cpu  versus  ()corr  at  all  instruments  for  (a)  good-SS  data  runs  ((5ss  <  0.1)  and  (b)  good- 
n(p)-R(p)  data  runs.  The  individual  data  points  are  represented  as  gray  dots  and  the  binned-means 
and  stds  are  shown  as  diamonds  and  vertical  bars,  respectively.  The  dashed  horizontal  line  at 
Cpu  =  0.9  indicates  the  Elgar  et  al.  (2005)  cutoff. 


Fig.  14.  (a)  and  (b)  R^'^^  versus  z'^^/as\g.  The  good-SS  data  runs  are  in  gray  and  the  good-/i^p^ 
data  runs  are  red.  In  (b),  the  dashed  horizontal  lines  indicate  the  cutoffs. 


Distance  from  the  Mean  Shoreline  (m) 

Fig.  1.  HB06  cross-shore  depth  transeet  versus  distanee  from  the  mean  shoreline  (blue  eurve). 
The  instrumented  frame  loeations  are  given  by  the  eireles  and  numbered  1-6.  An  additional 
instrumented  frame,  loeated  between  1-2,  was  often  buried,  and  is  not  eonsidered  here.  The 
typieal  tide  range  is  shown  with  the  horizontal  dashed  red  lines. 
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Fig.  2.  Example  ADV  measured  (a)  vertieal  veloeity  w,  (b)  signal-strength  (SS),  and  (e)  eorrela- 
tion  (CORR)  versus  time.  This  160  s-long  data  segment  is  from  instrument  1  (see  Fig.  1)  on  18 
September  0500.  The  red  dashed  horizontal  line  in  (b)  is  the  suggested  Elgar  et  al.  (2005)  eutoff. 
In  (e)  the  eorrelation  eutoffs  (  1)  for  /^  =  8  Hz  (7corr  =  0.526,  red-dashed)  and  the  mean  flow 
/s  =  0  Hz  (7CORR  =  0.3,  magenta)  are  shown.  The  water  depth  h  =  0.57  m,  Hsig  =  0.30  m, 
^adv  =  0.13  m,  and  z[-^  =  0.26  m. 


Fig.  3.  Fraction  of  bad  SS  data  points  6ss  versus  z[^/asig  at  instruments  1-6  (see  legend)  where  z[^ 
is  the  distanee  of  the  ADV  transdueer  below  the  mean  sea  surfaee  and  a^ig  is  the  signifieant  wave 
amplitude.  Note  that  no  instrument  6  data  points  are  present  in  this  axes  range.  The  horizontal 
dashed-line  is  the  (5ss  =  8  x  10“^  eutoff  for  disearding  a  data  run  (Elgar  et  al.  2005).  The  blaek- 
dashed  curve  is  the  proposed  scaling  (3)  based  upon  the  data. 
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Fig.  4.  Mode,  mean,  and  standard  deviation  of  bad  SS  gap  lengths  versus  5ss  at  all  instruments. 
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Fig.  5.  5corr  versus  5ss  at  all  instruments  (see  legend).  The  vertieal  dashed  line  indieates  the 
6ss  =  0.1  cutoff. 
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Fig.  6.  (a)  The  5corr  versus  ^^adv/^sig  instruments  1-5  where  log^o  of  the  energy-flux  gradient 
dF/dx  (m^  s“^)  is  eolored.  Only  good  SS  data  runs  that  satisfy  the  ^ss  <0.1  eriteria  are  shown. 
Data  runs  with  (5corr  >0.7  are  not  shown.  Panel  (b)  is  a  blow  up  of  the  nose  regions  in  panel 
(a). 
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Fig.  7.  versus  ^corr  for  good-SS  data  runs  at  all  instruments  (see  legend  in  Fig.  5). 


(a)  <5corr  =  0.028,  =  -1.67  ±  0.13  (b)  (5corr  =  0.255,  =  -2.14  ±  0.09 
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Fig.  8.  Example  of  vertical  velocity  spectra  Pww{f)  versus  frequency  /  at  instrument  3  that  are  (a) 
good  (time  =  240,  (5corr  =  0.028,  h  =  1.83  m,  Hsig  =  0.83  m,  ^^adv/^sig  =  3.32,  and  Cp„  =  0.86 
and  (b)  bad  (time  =  280,  5corr  =  0.255,  h  =  1.03  m,  Ffgig  =  0.57  m,  =  2.17,  and 

Cpu  =  0.91.  The  green  dashed  line  is  a  -5/3  power  slope,  (c)  and  (d)  are  the  e^"\f)  versus  /  over  a 
narrower  frequency  range  that  correspond  to  (a)  and  (b).  The  solid  black  line  is  e^‘^(/),  the  shaded 
region  are  the  error  bars  (derived  from  P/‘/)  and  the  blue  dashed  line  is  and  the  red  dotted  line 
is  the  linear  best  fit  slope. 


Fig.  9.  Vertical  velocity  spectra  P^w  power-law  exponent  (a)  and  (b)  versus  ^corr  for 
all  good  SS  data  runs  (5ss  <  0.1).  In  (c)  and  (d)  versus  (5corr  for  cases  that  additionally 
also  pass  the  —5/3-exponent  test  (7).  The  individual  data  points  are  represented  as  gray  dots  and 
the  binned- means  and  stds  are  shown  as  diamonds  and  vertical  bars,  respectively.  The  horizontal 
dashed  black  line  is  the  —5/3  slope.  Note  the  change  in  vertical  scale  between  (a,b)  and  (c,d). 
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Fig.  10.  Binned  number  N  of  good-SS  (iV,  eireles),  good-/i^p’  red  triangles),  and  good-fi^ 
iN^^\  green  diamonds)  data  runs  versus  (5corr- 


Fig.  11.  (a)  and  (b)  versus  (5corr  at  all  instruments  for  the  good  SS  (5ss  <  0.1)  data 
runs,  and  (e)  and  (d)  R^''>  versus  (5corr  for  good-fi'^'^^  and  good-/i^‘^  data  runs,  respeetively. 
Individual  data  points  are  represented  as  gray  dots  and  the  binned-means  and  stds  are  shown  as 
diamonds  and  vertieal  bars,  respeetively.  The  dash-dot  horizontal  line  represents  R  =  1.  The 
horizontal  dashed  lines  are  the  R  =  2  and  R  =  0.5  eutoffs,  respeetively.  Note  the  ehange  in 
vertical  scale  between  upper  and  lower  panels. 
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Fig.  12.  versus  5corr  at  all  instruments  for  (a)  good  SS  data  runs  (^ss  <  0.1)  and  (b) 

good-/i''’'-i?^'’’  data  runs.  The  individual  data  points  are  represented  as  gray  dots  and  the  binned- 
means  and  stds  are  shown  as  diamonds  and  vertieal  bars,  respeetively. 
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Fig.  13.  Cpu  versus  (5corr  at  all  instruments  for  (a)  good-SS  data  runs  (5ss  <  0.1)  and  (b)  good- 
^{p)-^{p)  (jata  runs.  The  individual  data  points  are  represented  as  gray  dots  and  the  binned-means 
and  stds  are  shown  as  diamonds  and  vertieal  bars,  respeetively.  The  dashed  horizontal  line  at 
Cpu  =  0.9  indieates  the  Elgar  et  al.  (2005)  eutoff. 


Fig.  14.  (a)  and  (b)  versus  -^adv/^sig-  The  good-SS  data  runs  are  in  gray  and  the  good-fi 
data  runs  are  red.  In  (b),  the  dashed  horizontal  lines  indieate  the  eutoffs. 


